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Possibility of Unconventional Pairing Due to Coulomb Interaction in Fe-Based 
Pnictide Superconductors: Perturbative Analysis of Multi-Band Hubbard Models 
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Possibility of unconventional pairing due to Coulomb interaction in iron-pnictide supercon- 
ductors is studied by applying a perturbative approach to realistic 2- and 5-band Hubbard 
models. The linearized Eliashberg equation is solved by expanding the effective pairing interac- 
tion perturbatively up to third order in the on-site Coulomb integrals. The numerical results for 
the 5-band model suggest that the eigenvalues of the Eliashberg equation are sufficiently large 
to explain the actual high Tc for realistic values of Coulomb interaction and the most probable 
pairing state is spin-singlet s-wave without any nodes just on the Fermi surfaces, although the 
superconducting order parameter changes its sign between the small Fermi pockets. On the 
other hand the 2-band model is quite insufficient to explain the actual high Tc. 
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Recent discovery of high-Tc superconductivity in iron 
pnictides has generated highly intensive research activi- 
ties in sohd state physics. After the discovery of supercon- 
ductivity in LaFeAsOi-a^F^: system,^) it has become evi- 
dent that transition temperature is raised above 40K by 
replacing La with other rare-earth elements (Ce, Pr, Nd, 
Sm, ...)^^^^ or by applying pressure.^-* Pairing mechanism 
and symmetry would be the most intriguing issues of this 
newly discovered superconductivity. In the present work, 
possibility of unconventional (i.e., not phonon-mediated) 
superconducting mechanism for Fe-pnictide superconduc- 
tors is investigated theoretically. There are already several 
reasons why unconventional pairing may be realized in 
iron pnictides: (i) Tc is high, compared with conventional 
(phonon-mediated) BCS superconductors, (ii) Fe3d states 
dominate the most part of the density of states near the 
Fermi level, and (iii) electron-phonon coupling is ex- 
pected to be weak by first-principles calculations. 

We introduce many-band Hubbard models for Fe3d- 
like orbitals and formulation. The Hamiltonian is given in 
the form H = Hq + H' . Hq is the non-interacting part: 
Ho = Eij E«' E<T Ui,ji'4g^Cjta. where cu^ (4^^) is the 
electron annihilation (creation) operator for Fe3d-like or- 
bital £ with spin a at site i. The tight-binding parameters 
tii,ji' are determined to reproduce a realistic electronic 
structure. H' is the on-site Coulomb interaction part con- 
taining four kinds of Coulomb integrals: U (intra-orbital 
repulsion), U' (inter-orbital repulsion), J (Hund's cou- 
pling), J' (inter-band pair- hopping). The same form of 
H' was used for the Ru4d£-like electrons of Sr2Ru04 in 
Ref. 12. Then the following linearized Eliashberg equation 
is solved numerically: 
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where a and a' are band indices, G'f\k) is the Green's 
function for band a (without self-energy corrections) , Ci 's 
are spin indices, VaaitT2,a'cT3a-4ik,k') is the effective pair- 
ing interaction, Aa,aa' (k) is the anomalous self-energy on 
band a, and A is eigenvalue. The effective pairing interac- 
tion is evaluated by third order perturbation expansion in 
H' . The third order perturbation theory has been applied 
to many other unconventional superconductors, and sug- 
gests correctly pairing symmetries, e.g., singlet d^2_y2- 
wave for cuprates and organic superconductors, triplet 
p-wave for Sr2Ru04, ... etc.^^' Transition point is deter- 
mined by A = 1. We take 32x32 k points and 512 Mat- 
subara frequencies for numerical calculations. 

Firstly, we adopt two-dimensional 5-band tight- 
binding model proposed by Kuroki et al.^"-' The electronic 
structure and the Fermi surface are given in Fig. 1 (in 
the unfolded representation, where each unit cell contains 
only one Fe atom). The Fermi surface consists of hole 
pockets around the (0,0) and (tt, tt) points and electron 
pockets around the (tt, 0) and (0, tt) points. In the orig- 
inal folded representation, where each unit cell contains 
two Fe atoms, the (0, 0) and (tt, tt) points are folded onto 
the r point, while the (tt, 0) and (0, tt) points onto the M 
point. The numerical results of eigenvalues are shown in 
Fig. 2(a). We see the most probable pairing symmetry is 
singlet s-wave and obtain sufficiently large eigenvalues to 
explain actual high T^s for realistic values of Coulomb in- 
tegrals {U = 1.2eV, U' = 0.9eV, J ^ J' = 0.15eV). is 
evaluated to be about lOOK. This is still higher than real 
values 20K-50K. If we include the self-energy corrections, 
then Tc will be decreased somewhat due to the effect of 
quasi-particle damping. The dx^-ya -electron component 
of anomalous Green's function Tx^-f^ (fe, iw„) is shown 
in Fig. 2(b) {X and Y axes are those in the original folded 
representation, while k = {kx,ky) is in the unfolded rep- 
resentation). The superconducting order parameter does 
not possess any nodes just on the Fermi surfaces, although 
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Fig. 1. Electronic band structure for the 5- band model. Inset 
shows the Fermi surface. 
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Fig. 2. (a) Eigenvalue A for various pairing symmetries as a func- 
tion of temperature in the 5-band model, (b) Anomalous Green's 
function of the local Fe3dj^2_y2-like orbitals, Fj(-2_y2 (fc, IttT), 
at T = 0.008. 



it changes its sign between the electron and hole pockets. 
In this sense the pairing symmetry is extended s-wave. 

We proceed to another model, i.e., 2-band model only 
for the dxz- and dy^-like orbitals, proposed by Raghu et 
al.") (See Fig. 3(a)). The maximum eigenvalue is given 
by the triplet p-wave pairing state, but is too small to 
explain the actual Tc, as we see in Fig. 3(b). Thus we 
may conclude that the 2-band model is quite insufhcient 
and the real multi-band situation is essential to describe 
the iron-pnictidc high- Tc superconductivity. 
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Fig. 3. (a) Electronic band structure for the 2-band model. In- 
set shows the Fermi surface, (b) Eigenvalue A for various pairing 
symmetries as a function of temperature in the 2-band model. 



In conclusion, our perturbation theory suggests that 
the iron-pnictide superconductivity may be unconven- 
tional one induced by electron correlation effect, as other 
unconventional superconductivity. The most probable 
pairing symmetry is s-wave without any nodes on the 
Fermi surface. One of the important differences from other 
unconventional superconductors is that the order param- 
eter will change its sign not on the Fermi surface but 
between the Fermi pockets. Possibility of triplet pairing 
is excluded. 
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